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Abstract 


Normal mode theory, extended to the slightly laterally heterogeneous Earth 
by the flrst-order Born approximation, is applied to the waveform inversion of 
mantle Love waves (200-500 sec) for the Earth’s lateral heterogeneity at 1=2 and 
a spherically symmetric anelasticity (Q^) structure. The data are from the Glo- 
bal Digital Seismograph Network (GDSN). The 1=2 pattern is very similar to the 
results of other studies that used either different methods, such as phase velo- 
city measurements and multiplet location measurements, or a different data 
set, such as marfitle Rayleigh waves from different instruments. The results are 
carefully analyzed for variance reduction and are most naturally explained by 
heterogeneity in the upper 420 km. Because of the poor resolution of the data 
set for the deep interior, however, a fairly large heterogeneity in the transition 
zones, of the order of up to 3.5% in shear wave velocity, is allow^ed. It is 
noteworthy that Love waves of this period range can not constrain the structure 
below 420 km and thus any model presented by similar studies below this depth 
are likely to be constrained by Rayleigh waves (spheroidal modes) only. 

The calculated modal Q values for the obtained model fall within the 
error bars of the observations. The result demonstrates the discrepancy of Ray- 
leigh wave Q and Love wave Q and indicates that care must be taken when both 
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Rayleigh and Love wave data, including amplitude information, are inverted 
simultaneously, 

Anom,alous amplitude inversions of G2 and G3, for example, are observed for 
some source-receiver pairs, This is due to multipathlng effects. One example 
near the epicentral region, which is modelled by the obtained I =8 hetero- 
geneity, is shown. 


1. Introduction 


Putting constraints on the lateral heterogeneity of the Earth is one of the 
most important contributions seismologists can make to the understanding of 
the dynamical behavior of the Earth. Since the work of Toks’dz and Anderson 
(1966), most studies of the global scale heterogeneity of the upper mantle have 
used surface waves. The most recent contributions using this approach are by 
Nakanishi and Anderson ( 1983a, 1983b). Recently, Silver and Jordan (1981) and 
Masters et al, (1982) used standing wave analysis (free oscillations), which is, in 
essence, not so different from phase velocity measurements, In an entirely 
different approach, Woodhouse (1983) inverted the observed seismograms 
directly for the heterogeneity of the Earth. His result showed a pattern similar 
to that found by Masters et al, (1982). His derivations of the formula used for 
the inversion were quite complicated, however. Tanimoto (1983) showed that an 
application of the first order Born approximation to the conventional normal 
mode approach (Gilbert, 1970) gives equivalent formula. Recently, with a 
slightly simplified approreb, Woodhouse and Dziewonski (1983) obtained a 
heterogeneous Earth model with the maximum angular order f=8, using prob- 
ably the largest quantity of data so far. 

In this paper, we apply the Born approximation method to mantle Love 
waves (200-500 sec) obtained from GDSN (Global Digital Seismographic Network) 
stations. We invert the seismograms directly for the lateral distribution of rigi- 
dity for different depths. The results are similar to those obtained by Masters et 
al. (1982), Woodhouse (1983), and Nakanishi and Anderson (1983a, 1983b). Con- 
sidering the differences in the methods and/or the data, the similarity is quite 
remarkable. In the process of inversion, we also perturb the cpherically sym- 
metric anelasticity parameters Q^, because no recent Earth model adequately 
explains Love wave Q values. Most recent Q models are constructed for Rayleigh 
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wave Q^ and there seems to be a yet unexplained difference between Rayleigh 
wave Q and Love wave Q at present. 

We also show an interesting multipathing effect near the source region, 
which can be explained by the heterogeneity obtained. 


2. Method 

The equation of motion is given by 

(po + <5p) a?u + (Ho + H)u =/(t,X 3 ). (1) 

where the subscript 0 refers to the zeroth order spherically symmetric Earth, 
6p and H are the deviations from it, and/(t,Xg) is a source as a function of time t 
and position x,. Tanimoto (19B3) showed that, under the first order Born 
approximation, the seismograms in a slightly heterogeneous Earth can be writ- 
» ten as 

u = Re 2 ^ + E 

k '^k , tn m'm -'^k 


+ 


3.ki»s ^k 


m'm 


FJ"u]^°^'”’(<lcm'| H|sm>— (u|<km'[(5|PS sm>) 


+Fi^'ui°>'"(<sm 1 H 1 km'>-«^<sm 1 6p | km'>)V 


( 2 ) 


where is an eigenfunction of a k-th multiplet in the zeroth-order spheri- 

cally symmetric Earth and satisfies 


(4 4°)*^ = Ho 


<km| 6p I sm'> = /I <5p 4°**” 4°^*^' 

•/ g 
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<kmjH|sm'> s aV, 

Ff = /j. ni")“7{p,x,) dV, 


and 


«k = + 


2 uil°)'"'(<km' ! H ! km>-w|<km' | Sp | km>)Fi?' 
m'm 

m 


Here 7(PiXg) is the Laplace transformation of /(t.Xg) with respect to time. This is 
equivalent to the formula derived by Woodhouse (1983) in a much more compli- 
cated way. The first term with the sum over m in [] represents a seismogram in 
a spherically symmetric Earth, and the rest are due to the lateral heterogeneity 
of the Earth. The second term in [] with the sum over m* and m is the contribu- 
tion from the same multiplet k, and the rest, i.e. all terms in the second and the 
third lines, are the contribution from different multiplets s (5«ik). Generally 
speaking, <km'|H|sm> and <km'|(5p|sm> are much smaller than <km'jH|km> 
and <km'|(5p|km>. Thus the last terms are small unless CJki=3Wg, that is when 
there exist modes with close eigenfrequencies. This condition is satisfied for the 
frequency range considered except for a few modes, For this reason, we drop 
the last terms in the following analysis, although it should be interesting to 
study these exceptional cases. 

In the formal Born approximation approach, the second term is not exactly 
a contribution from the same inultiplat k, but a contribution from modes with 
close eigenfrequencies. The last terms are contributions from modes with dis- 
tant eigenfrequencies. If coupling of modes occur, (Luh, 1973, 1974; Woodhouse, 
1980; Morris and Geller, 1982; Tanimoto and Bolt, 1983; Masters ei al., 1983), ,it 
should be included in the second term. Summation must be extended to other 
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multiplets in addition to the k-th one in such a case. Formally, the applicability 
of Born approximation approach is controlled by the magnitude (smallness) of 
terms associated with H and 6p compared with terms in the zeroth-ordef spheri- 
cally symmetric Earth. However, in the following applications in this paper, the 
coupling effects are ignored. 

In the geometrical optics approximation, Jordan (1978) showed that the fre- 
quency shift of the multiplet is the average of the local perturbation to the 
eigenfrequency over the great circular path. Similar proof can be established for 
attenuation in the same manner, i.e. the Q of the multiplet is the average of the 
local perturbation to the attenuation over the great circular path. Since this is 
not related to other topics in this paper, it is shown in the Appendix. 

We find it convenient to change from, epicentral coordinate, as In (2), to geo- 
graphic coordinate, as Woodhouse and Girnius (1982) did. Then after dropping 
the terms in the second and the third lines, (2) becomes 

u = Re 2 ^ ( 2Rk"Sr - -Vs , (3) 

k "C m ui Xi 


where 


N=-l 


= - £ Fj? 

N =-2 


and 


Y^^(e,<fi) = PjN'^(cos^)e^'^® 

is the generalized spherical harmonics of Phinney and Burridge (1973). Here 6 
and <p are colatitude and longitude in the spherical coordinate and the 
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expressions for F|^ are given in Woodhouse and Glrnius (1982), 

Computation of is done by using the recursive reiation for 

which is 

= -a„a pf^ P/^ 

with 

_ _ /TiTm+^rai-ry 

V (f4-m+l)(f-m) 

and 


The recursion is started from, m=i, noting 

Pl^Kcos^ =.' (-1)‘-N>y (cos|-)'-"N (sin|-)‘"N 

This can be found in Edmonds (1960), whose notation is instead of 

Pi^^(cosR). 

Compiete expressions for HiF*" are given in Woodhouse and Dahien (1978) in 
the absence of coupling effects, In our particular case, this becomes 

3_2 ]' 

fi + 01Pn''^k i — J^yjTjy *^in'm 

a 

^ f Ms r^dr, (4) 

S 0 

where Q is the sidereal rotation rate, 7.292115 x 10~® (1/ s), ae^ is the ellipticity 
splitting parameter for a given multiplet k (Dahien, 1976; Woodhouse and Dahien. 
1978), 
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a the radius of the Earth, the spherical harmonic expansion of aspherl- 

cally distributed rigidity (the angular order S, the azimuthal order m'-m), and 


2wM, = 


K^+1-) - |'s(s+l) }{ 


dvr w vg 
8r r ' 


■h 


1 ( 1 + 1 ) 


1 ( 1 + 1 ) - 2 


|-s(s+l) 


41(1+1) - s(s + l) - 2 



where 


s l"\ 

Q Q q] is the Wigner’s 3j symbol (e.g. Edmonds, 1960) and w is an eigen^ 


function of a torsional mode. We used the eigenfunctions of 1066A in this study 
(Gilbert and Dziewonski, 1976), 


The first term in (4) is the contribution from rotation, the second from the 
hydrostatic ellipticity, and the third from aspherical distribution of rigidity. We 
dropped the density perturbation and the boundary perturbations in this study. 
We subtract the synthetic seism.ogrv\ms for the reference spherically symmetric 
Earth (1066A; with the rotation and ellipticity effects from the data and invert 
the resultant disturbances for 

We found, however, that amplitudes were not fitted very well if we use the 
recently published Q models. This is not surprising, since most Q models are 
constructed in order to satisfy Q of Rayleigh waves and there seems to exist, at 
least at present, an inconsistency between Love wave Q and Rayleigh wave Q. 
Thus we adopted the Q model of Masters et al. (1983) as the starting model and 
perturbed them. Partial derivatives of the seismograms with respect to the 
intrinsic are easily obtained, since Q“^ of surface waves are linear functionals 
of (Anderson and Archambeau, 1964), 
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In the radial direction, we use two shells, one between the surface and 420 
km and the other between 420 km and 670 km. In each shell, we solve for 
(5/Uol<5/i|',5/z^^^(5/x^^and<5qm(• 5(1/ Qp)). Thus In total we have 14 parameters. 
Conventional least-squares Is used In solving this. 


3. Data 

We analyzed Love waves from 11 large earthquakes In 1980 recorded by 
GDSN stations. A list of earthquakes and the source parameters used are given 
in Tables la and lb, We used the source parameters determined by Nakanishi 
and Kanamori (1983) and did not perturb them during the inversion, since Love 
waves can place only four independent constraints on moment tensors. 

In order to reduce the size of the data set, the seismograms were filtered 
between 2 mHz (500 sec) and 5 mHz (200 sec) and resampled at 40 sec intervals, 
A cosine taper was applied between 2 and 3 mHz and also between 4 and 5 mHz. 
An example of the data set for the earthquake in northern California (No. 24 In 
Table la) is shown in Fig. 1. G1 is not used in the analysis, and if G2 and G3 were 
clearly recorded the data were analyzed further, Fig, 1 shows some of the best 
examples in our data set. Tn many seismograms, glitches were present and the 
seismograms had to be truncated, In total, we found 60 source-receiver pairs 
with excellent S/N ratios. The great circle paths are showm in Fig. 2, 

4. Analysis 


First, we calculate partial derivatives of seismograms with respect to 
5q(= <5(1/ Q)) and <5/zf* for each source-receiver pair. Equation (3) can be writ- 
ten symbolically as 


auo . xn ^ 9^0 




( 5 ) 


where the subscript I refers to the shell number In the radml direction and <J/U|7 
is the (l,m) component of spherical harmonic expansion of 

Instead of solving for complex w!:iu;h Is the case when we use complex 
spherical harmonics we changed to P/" (cos <9) cosmos and 

Pt'” (cos(9) sinm^ at this stage. Formulation up to this point is done by Yf* {6,p) 
of Edmonds (1960), It is advantageous to use them in the formulation, since we 
can fully use the properties of 3j symbols, but in computation real quantities are 
much easier to handle. UVe normalize the associated Legendre polynomial by 
[(2-dmo)(2i + l)(f-nri)l/(Z+m)!j''^^ and use the notation Pp (cos0). 

Fig. 3a shows an example of a set of partial derivatives for one source- 
receiver pair. The earthquake in Fig. 3 is No. 1 in Table la, and the receiver is 
ZOBO in Bolivia. All fundamental, first overtone, and second overtone modes 
between the periods of 200 sec and 500 sec are included in these calculations. 
This is the most time-consuming step in the inversion, taking about 20 hours on 
the VAX 11/780 computer. 

In Fig. 3a, we show the data, the synthetic seismogram with ellipticity and 
rotation corrections, their differences, partial derivative with respect to 1/ Q in 
the upper 420 km and that in the transition zone and partial ucrivatives with 
respect to corresponding spherical harmonic components of rigidity in the 
upper 420 km and those in the transition zone. Four large amplitude wave pack- 
ets are G2, G3, G4 and G5 from left to right. The small amplitude oscillations in 
between are contributions from overtones. Contributions of overtones are 
larger in the deeper transition zone than in the upper 420 km, but in both shells, 
the fundamental modes are dominant. Fig. 3b is the case of a different source- 
receiver pair, which shows the same features. From these two figures, we expect 
that results are mostly determined by fundamental modes. 
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The form of the problem becomes 

A X = b (6) 

where the elements of A are given by -^^and 

3qj ddfMf 

= {6qi,6q%,6fi§, * • • ) 

and the elements of b are given by u-Uq at each sampled point In the selsmo- 

4 

grams. Simple least stiuares method is used to estimate x, i.e. 

X = {^JA)‘*‘A'^b . (7) 

Standard error of i-th parameter in x, ci, is estimated by 

"• = V £ ( S - b,)= . (8) 

V n-m ^ 

where (A'^A)"Hl is the i-th diagonal element of (A'^A)"S n the number of data 
points, m the number of parameters used for fitting, and an^, x^, and b^ are ele- 
ments of A, X, and b, 

5, Results 

Heterogeneity 

We performed inversions for five cases as shown in Table 2. The cases differ 
in whether Q is included and whether rigidity in the upper region (0-420 km) 
and/or in the lower region (420-870 km) are included. In the table, o means 
that the parameter is included, while x denotes that it is not. When Q is 
included, in both the upper and the Lower regions are perturbed. We calcu- 
lated the variances for the initial model (spherically symmetric Earth with 
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rotation and elliptlcity correction), erf, and for the resultant model, cr|. Variance 
reductions are (of - af/ of and arc tabulated in Table 2. 

Two features of this table deserve particular attention, First, as comparis- 
ons of 1(42.6%) with 2(42,2%) and 3(21,7%) with 4(21.2%) make clear, the Interac- 
tion of the perturbation in the lower region achieves only an extremely small 
improvement in variance reduction. Second, Incorporating Q In the inversion 
produces large variance reductions, This point is discussed in the next section. 

Inverted results for the upper region are shown 'in Fig, 4 for each case, 
« 

Coeflficients of 1=2, used for plotting, are given in Tables 3a and 3b, The chain 
and dashed lines correspond to positive and negative regions, respectively, and 
the solid line is the zero line. The contour interval is 0.4% in this figure. The 
similarity of the four figures suggests that the results of the aspherical distribu- 
tion of 6/1 in the upper region are almost independent of the introduction of Q 
and also the lower region in the inversion. Only slight shifts of the patterns are 
observed. The coefficients of f=2 are tabulated in Table 3 with the estimated 
standard error for each case of inversion. 

It is interesting to note, however, that case 5, in which all lateral hetero- 
geneity is confined to the transition zone, gives smaller but comparable variance 
reductions (19.1%). The results are shown In Fig. 5, along with the results of 
cases 1 and 3 for the lower region, The results of case 5 display a very similar 
pattern to the results in the upper mantle shown in Fig. 4, but contour intervals 
are 5 times larger (2%), The peak value, 6 /i Is about 7 %, as compared to about 
1.5 % in Fig, 4, This can be understood from Figs 3a and 3b. Partial derivatives 
for the lower region have virtually the same patterns as those for the upper 
region, but the amplitudes are about 5 times smaller. This is also the reason 
that the extra variance reductions achieved for cases 1 and 3 are so small com- 
pared to cases 2 and 4. Since the shapes of the partial derivatives are similar, 


- 10 - 


Ihe extra degrees of freedom Introduced into the model space by the introduc- 
tion of parameters in the lower region are very small, And since the amplitudes 
are 5 times smaller, the extra variance reductions are very small. 

This can possibly be avoided by using deep earthquakes, because as shown 
in Figs 3a and 3b the amplitudes of partial derivatives due to overtones tend to 
be greater for the lower region than for the upper region, which is opposite to 
the trend of fundamental modes. Thus, for deep earthquakes, the shapes of the 
partial derivatives can be dliTerent for different depths. The earthquakes in this 
paper, unfortunately, are all from shallow depths. Ifaturally, the resolution of 
our data set in the lower region Is not good and our data set does not require 
lateral heterogeneity In the transition zone, But at the same time It can not rule 
out the pn&'^lbiUty of a heterogeneous transition zone (420-670 km). 

6 compares the results in the upper region in case I with the results of 
Nakanishi and Anderson (1983) and Woodhouse (1983). The results of Nakanishi 
and Anderson are for Love wave phase velocities at the period of 307 seconds, 
which corresponds approximately to the dominant periods in our data set, The 
results of Woodhouse (1983) are obtained from Rayleigh waves using IDA (Inter- 
national Deployment of Accelerographs) network. Considering the differences in 
methods (waveform matching and phase velocity measurements) and in data 
(Love waves from GDSN and Rayleigh waves from IDA), the similarity of the 
results is quite remarkable. The results of Masters et al, (1982), which uses yet 
another approach (multiplet location measurements), also show a very similar 
pattern. We should note, however, that the interpretation of the depth of 
heterogeneity is different among these studies. Masters et af,(l982), for exam- 
ple, interpret their results entirely in terms of heterogeneity in the transition 
region. 
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Attonualioa 

The results for are discussed here. In Table 4, the initial model (Masters 
et al., 1983) and the results of this study are shown. Differences between them 
are quite large, which confirms nothing but the inconsistency of Rayleigh wave 
and Love wave Q Values. Most recent Q models are constructed mainly from Ray- 
leigh wave (spheroidal mode) data. This is important because of its implication 
that simultaneous waveform Inversion of Rayleigh and Love wave data has to be 
done with care. A Q model that fits Rayleigh wave data will give systematically 
higher amplitudes for Love waves and vice versa. 

Q of torsional modes for case 1 are calculated and shown in Fig. 7 with 
observed torsional Q values reported in a table of PREM (Dzlewonski and Ander- 
son, 1902), The error bars are for observed Q valU'^’i. as reported there. The 
symbols o and x correspond to calculated fundamental and first overtone 
modes. Q values of fundamental modes are rather low, while those of first over- 
tones are not. Overall, Q values obtained in this study are consistent with 
reported Q values in PREM. 

There is one caution, however. The source parameters are determined by 
Rayleigh waves using a certain Q structure (Nakanishi and Kanamorl, 1983). 
Estimated moments depend upon the Q model, Since those moments are used • 
in the process of inversion, there can be some trade-off between them. The fact 
that 60 paths for 11 earthquakes are u.sed may have helped to reduce the trade- 
off, but there can still remain a systematic bias in estimated Q values. 



6. Discussioa 

Figs. 8, 9, and 10 show the comparisons of waveforms before and after the 
inversion in case 1. In each figure, the top seismogram Is the data; the second 
the synthetic seismogram for 1066A with the Q model of Masters p.t al. (1903) 
and the rotation and ellipticity correction; the third the same with the top one 
(the data); and the fourth the synthetic seismogram for our laterally hetero- 
geneous Earth model. 

Fig. 8 shows an interesting r^iultipathing effect. The earthquake is No. 24 in 
Table ia and the receiver is in Albuquerque, New Mexico (U.S.A,), Since the 
receiver is close to the source (about 15®), G2 and G3 and also G4 and G5 arrive 
at the receiver closely in time. The interesting feature in the data is the larger 
amplitude of G3 than that of G2, which is not modelled in the spherically sym- 
metric Earth. However, this feature is explained in the seismograms after the 
inversion. Thus the obtained f=2 heterogeneity is capable of producing this 
multipathing effect, although the heterogeneity is not so strong. This is because 
the receiver Is close to the source, where the effect of heterogeneity appears 
strongly. This should also be true for receivers near anti-podes, but there w'ere 
no anti-podal data in our data set. 

Fig. 9 shows the seismogram at Albuquerque for earthquake No. 8. This 
seismogram also has G2 and G3 in the first wave packet and G4 and G5 in the 
second wave packet. In this case, the main improvement is due to anelasticity. 

Fig. 10 shows the worst example in the inversion. The source is No. 4 in the 
Kuril Islands, and the receiver is in India. The fit of the seismogram clearly 
became worse after the inversion. In particular, the arrival of the wave packet 
is too slow and is apparently perturbed in the wrong direction by the hetero- 
geneous model. In order to see the point, the time interval from left to right is 
taken to be 4 times as large as that in .Figs. 8 and 9. This is presumably because 
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i-2 heterogeneity is not enough to represent the Earth’s heterogeneity. L=2 
heterogeneity should he large, since fairly large variance reductions are 
achieved and the fit of -the synthetic seismograms to the data is improved in 
most seismograms. However, an example like Fig. 10 shows that spherical com- 
ponents of lateral heterogeneity for should counteract 1=2 heterogeneUy 
for some source-receiver pairs. 

Although we do not discuss this problem extensively, there may be a prob- 
lem of aliasing. This problem arises because of the existence of short wavelength 
heterogeneity in the Earth. When the maximum I In the inversion is not high 
enough, these short wavelength oscillations can affect the estimation of low 
order coefficients. There does not seem to be a simple cure for this problem, as 
long as spherical harmonic approach is used. We will address this problem in a 
later contribution. 

7. Conclusion 

Using mantle Love waves from the GDSN network, we obtained the lateral 
heterogeneity of f=2 pattern and spherically symmetric by waveform, inver- 
sion. Partial derivatives are calculated by using the first Born approximation to 
the normal mode approach. The obtained 1=2 pattern is very similar to the 
results of Masters et al. (1982; multiplet location measurements), Woodhouse 
(1983: waveform-inversion of Rayleigh waves from the IDA network), and Xakan- 
ishi and Anderson (1983; phase velocity measurements). Our results can most 
naturally be explained by heterogeneity in the upper 420 km. Because of the 
poor resolution in -the transition zone, the data allow fairly large heterogeneity, 
of the order of 3.5% in shear wave velocity, in this region, 

structure is simultaneously inverted, and a reasonable model, which 
explains the observed Q of torsional modes, is obtained. Howe> , possibility 
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of bias due to the trade-off between the Q structure and the estimated moments 
of the earthquakes remains, It is important to include Q in the nversion since no 
recent Q model satisfy torsional modes and a fairly large additional variance 
reduction is achieved by incorporating them. This result conflrmr the incon- 
sistency of Rayleigh wave Q and Love wave Q values. Thus care must be taken in 
performing waveform inversion of Rayleigh waves and Love waves simultane- 
ously, 

Observed muitipathing near the source is quite well modelled by the 1-2, 
heterogeneity obtained in this study. But there are a few examples in which the 
fit of the synthetic seismograms becomes worse. This simply means that the 
Earth’s heterogeneity is more complicated than i=2 pattern obtained. 
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Appendix: Q of a muiUpict Is the average ot the local perturbation to the 
attenuation over the great circular path 

The relation between the elgenfrequency of a mulLip.et in a heterogeneous 
Earth, 2Jk< and that in a zeroth-order spherically symmetric Earth, Is given by 


2 aii°^'”'(<km'|H!km> - w^<km'!dp!km>)Ff 
Wk = Wk + (A-1) 

2wk E 4°^"’ 

m 

Attenuation can be incorporated in this formula by treating the elastic con- 
stants K and in H as complex. We write Wk as 

Wk(0.'^) = Wk +■ \k(0.‘^) -t iAk(0,^), (A-S) 

where (0,'l' is the pole of the great circular path under consideration. 

Jordan(l978) introduced the quantity, dwiocgj(0,$), which is the local pertur- 
bation to the eigenfrequency, and defined it by 

a 

<5«ioca](®*^) = /M*^(r)dm(r,0,5i>)r®dr, (A-3) 

0 

where <5m is the perturbation of elastic constants and the corresponding ker- 
nels. Then he showed 


Xk(e,$)^ E Ps(0) £ (<5W!ocal)s^Y3‘0.$) (a-4) 

3=S, even t=-s 


where 


<5W!ocal(^.^^’) = E E(‘5‘^local)s YK^,»5). 

. 9 t 


For attenuation, we define 
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5u\d,<p) = / 
0 


K^(r) 6ic{Ti6,(p) + U'^(r) 


dr 


(A-5) 


where Sk and d/j, are the imaginary parts of k and fx, and are corresponding 
kernels. Note that the local Q value, QiocaiC^i^’). is given by 

1 _ 26uH9,(fi) 

Qlocal(^»i^) *^k 


In exactly the same way as Jordan, we can obtain 

f 

Ak(S,#) « S P,(0) t (5o'')>Y‘(a,$). (A-6) 

gs3 t=-a 


where 


6oHe.v) = E S 

s t 


Using the relation (Backus, 1964, equation 44), 




and extending the upper limit of the summation, 21, to infinity (geometrical 
optics approximaton; Jordan, 1978), we obtain 


Ak(9.‘^) = 2^ A 


1 


A 


ds 


2 2rra -^l 


(A-7) 


where ds is the line integral along the great circular path. Q of a multiplet is 
related to Ak by 


Wk 


2Qa^3(e.<l) ‘ 


Ak(0,$) = 


thus we have 
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1 1 r ds 


(A-B) 


Note that this is the relation between the standing wave Q (temporal Q) and 
for propagating wave Q (spatial Q), a factor of (group velocity)/(phase velocity) 
must be multiplied as weights. 
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FIGURE CAPTIONS 

Fig. iJ An example of filtered seismograms {2 - 6 mHz) for the earthquake No. 24 
in Table la, 

Fig. 2; Great eircle paths of 80 source-receiver pairs. 

Fig. 3a: An example of a set of the data and the partial derivatives for one 
source-receiver pair. The earthquake is No. 1 in Table la and the receiver 
Is ZOBO in Bolivia. The top three traces are the observed data, the syn- 
thetic seismogram for the model 1066A with rotation and elliptlclty 
correction and the residual between the two. The next two are the partial 
derivatives with respect to q (= 1/Q) in the upper 4-20 k,m and in the tran- 
sition zone. The next six traces are the partial derlvativos with respect to 
the spherical harmonic components of 5//, in the upper 420 km. For exam- 
ple, (21C) denotes the coefficient of P 2 ‘(cos^ cos^. The last ?ix traces are 
those for the transition zone. 

Fig. 3b: Same as Fig. 3a except for a different source-receiver pair. 

Fig. 4; The results of the inversion in the upper 420 km. The numbers on the top 
right corner correspond to the case in Table 2. L=2 components of are 
plotted at an contour interval of 0,4 %. The dashed lines correspond to 
negative region, the chain lines to positive region and the solid lines are 
the zero lines. The negative peak is about 1.5 % both on the East Pacific 
Rise and India. In terms of shear wave velocity, this becomes about 0.7%, 

Fig. 5: The results of the inve^'sions for the transition zone. Note that the contour 
intervals are larger than those in Fig, 4. 

Fig. 6: Comparison of the results of case 1 with the results of other studies. The 
second is Love wave phase velocity taalysis by Nakanishi and Anderson 
(1983b) for the period of 307 sec. The third is the. result of Wood- 
house(l903), using a similar method but Rayleigh waves from IDA net- 
work. 

Fig. 7; Computed torsional mode Q values for the inverted model. The sym- 
bols o and X correspond to fundamental and first overtone modes. The 
symbols • with error bars are the observed torsional Q values found in 
Dziewonski and Anderson(l902). 

Pig. 8: Comparison of the data and the synthetic seismograms before and after 
the inversion. The first and the third traces from top are the observed 
seismogram, the second the synthetic seismogram for 1066A with rota- 
tion and ellipticity correction, the fourth the the synthetic seismogram 


I 
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for the obtained heterogeneous Earth. In this example, the receiver is 
cicse to the source (about 16 degrees), so the first wave packet contains 
both G2 and G3 and the second G4 and G5. Amplitude inversion of G2 and 
G3 is reproduced in the synthetic for the obtained heterogeneous model. 
Pig, 9: Same as Fig.B except for a dlfterent source receiver pair. 

Fig. 10: The worst example In our inversion. Clearly, the wave packet Is per- 
turbed in the wrong direction by the obtained heterogeneity. This 
demonstrates a rather trivial point that the Earth's heterogeneity can not 
be represented accurately by 1=2 term only. 


Table la 

I4st of Earthquakes (in 1900) Used 


No. 

m 

d 

h 

m 

s 

Lai. 

Lon. 

Ms 

Refflon 

1 

1 

1 

16 

42 

40.0 

30.01SN 

27.7B0W 

6.7 

Azores 

3 

2 

7 

10 

49 

16.0 

54.150S 

15B.09OE 

6.5 

Macquarie Islands 

4 

2 

23 

5 

51 

3.2 

43.530N 

146.753E 

7.0 

Kurile Islands 

7 

3 

24 

3 

59 

51.3 

52.969N 

167.670W 

6.9 

Fox Islands 

8 

6 

9 

3 

28 

10.9 

32.220N 

114.9B5W 

6.4 

Cal-Mex Border 

16 

7 

29 

3 

11 

56.3 

13. 10 IS 

166.338E 

8.7 

Vanuatu Islands 

19 

10 

10 

12 

25 

23.5 

36.195N 

1.354E 

7.3 

Algeria 

24 

11 

8 

10 

27 

34,0 

41.117N 

124,25317 

7.2 

N. California 

26 

11 

23 

16 

34 

53,0 

40.914N 

15.368E 

6.9 

Italy 

27 

12 

17 

16 

21 

58.8 

49.479N 

129.498W 

6.8 

Vancouver Island 

28 

12 

31 

10 

32 

11.0 

46.060N 

151.453E 

6.5 

Kurile Islands 
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Table lb 

Source parameters of the earthquakes used 


No. 

depth(km) 

T 

Mo(xlO^O) 

5 

X 


1 

9.75 

17.2 

2.38 

86.2 

3.0 

-31.0 

3 

9.75 

29.7 

1.9 

84.0 

0.0 

-70.0 

4 

43.0 

19,3 

6.31 

70.0 

39.2 

27.0 

7 

33.0 

30.1 

2.95 

60.0 

80.1 

53.3 

3 

9.75 

15,4 

0.465 

90.0 

180.0 

140,1 

16 

43.0 

19.1 

1.57 

54.0 

93.5 

160.0 

19 

9.75 

30.2 

4.89 

54.0 

81.8 

225.0 

24 

16.0 

31.7 

10.3 

'90.0 

0.0 

49,8 

26 

9,76 

44.7 

2.84 

63.0 

275,0 

-43.0 

27 

9.75 

26.2 

1.54 

90.0 

180.0 

.37 1 

w 1 • .a. 

28 

33.0 

27.8 

2.90 

68.0 

89,6 

28.3 


KEY; 

T - rise time \ec) 

Mo " seismic moment (dyne cm) 

6 - dip angle (deg) 

X - slip angle (deg) 

- strike, measured clockwise from north (deg) 
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Table 2 

Results of Inversion for 5 cases with variance reductions 



reduction (%) 42.6 42.2 


21,7 


21.2 


19.1 
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Table 3a 

Results of the inversion for f =2 coefficients of 
Sigma (a) is the standard error calculated from (0). 


Upper 


Lower 



case 1 


case 2 




(X ■ 

W' 

a 

20 

O.llB 

0.060 

0.045 

0.017 

21C 

0. 152 

0.027 

0.107 

0,007 

21S 

-0.093 

f 

0.058 

-0.196 

0.015 

22G 

0.606 

0.043 

0.615 

0.011 

22S 

-0.073 

0.056 

-0.353 

0.014 

20 

-0.397 

0.353 



21C 

-0.239 

0.141 



21S 

-0.563 

0.305 



22C 

0.045 

0.223 



22S 

<1.495 

0,291 
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Upper 


Lower 


Table 3b 



case 3 


case 4 



tf/ur 

a 


or 

20 

-0.042 

0.079 

-0.079 

0.020 

21C 

0.114 

0.032 

0.108 

0.008 

21C 

-0.068 

0.068 

-0,191 

0.017 

22C 

0,639 

0.050 

0.616 

0.012 

22S 

-0.049 

0.085 

-0.351 

0.017 

20 

-0.203 

0,411 



21C 

-0.037 

0.165 



21S 

-0.667 

0.356 



22C 

-0.110 

0.260 



22S 

-1.812 

0.340 




case 5 






a 



20 

-0.336 

0,106 



21C 

0.525 

0.042 



21S 

-0.977 

0.090 



22C 

3.104 

0.066 



22S 

-1.B32 

0.087 




Lower 
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Table 4 

Starting model(Masters et al., 1903) and the inverted results 



Masters et al.(l983) 

Case 1 

Case 2 

0 - 430 km 

120 

91 i 1 

91 i 1 

420 - 670 km 

280 

168 ± 20 

163 ± 20 

670 - 2887 km 

340 

Not perturbed 

Not perturbed 
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